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\ Abstract. Lee, in a series of papers, described a unified formulation 

CN ■ of the statistical thermodynamics of ideal quantum gases in terms of 

the polylogarithm functions, Li s (z). It is aimed here to investigate 
| the functions Li s (z), for s = 0, —1, —2, . . . , which are, following Lee, 

referred to as the polypseudologarithms (or polypseudologs) of order 
n. Various known results regarding polypseudologs, mainly obtained 
in widely differing contexts and currently scattered throughout the 
literature, have been brought together along with many new results 
and insights and they all have been proved in a simple and unified 
manner. In addition, a new general explicit closed-form formula for 
these functions involving the Carlitz-Scoville higher tangent numbers 
has been established. 
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1. Introduction 

Lee, in a series of papers [THS], not only has shown for the first time that statistical 
thermodynamics of the Bose-Einstein and Fermi-Dirac ideal gases can be unified, 
but also has obtained a unified description of the main thermodynamical functions 
in terms of the polylogarithmic functions Li s (z) (for their definitions see Sec. 2) and 
also examined details of certain physical phenomena revealed by this approach. 

The unification is based on the following expression for the reduced density of 
non-relativistic ideal gas of N particles confined to a ti-dimensional box of volume V 

pX d tj.\-r- />\ j- f z if Bose-Einstein , . 

— = sgn(C) Li, /2 (C), C =(_, if Fermi-Dirac, (L1) 

where p = N/V = N/L d is the particle number density, A = (27r/3/m) 1 / 2 (with 
h = 1) is the thermal wave-length, (3 = l/ksT and g = 2s + 1, while &b, T, m 
and s denote the Boltzman constant, the temperature, the mass and the spin of the 
particles, respectively. Finally, the fugacity z is related to the chemical potential /i 
of the system as z = e^ M . The grand partition function then easily follows from (1.1) 
Q = exp {sgn(C)(A/L) d Li rf / 2 +i(C)} > while the basic thermodynamic functions, such 
as the energy U, the entropy S, the pressure P and also fluctuation quantities such 
as the number fluctuations Y, can be expressed in polylogarithms because of their 
relationship with the density p. 
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The important property of the unified formulation is that U, S, P, Y and the 
reduced density depend only on two variables, the fugacity and dimensionality, and, 
more importantly, their values are given by the very same function - polylogarithms 
(of order s and argument z), Li s (z); the fugacity z determines the argument, while 
the dimensionality d enters the polylogarithms through its order-an integral order if 
d is even, and a half-integral order if d is odd. Moreover, by the principle of analytic 
continuation, the obtained results may be assumed valid for any d. 

It is noteworthy that the polylogarithm formulation of statistical thermodynamics 
evidently reveals that the thermodynamic properties of ideal quantum gases are 
described by the structural properties of Li s (z) and, therefore, there is a need to 
understand better these important but little known functions. In addition, it is 
not surprising that this formulation is a quite natural framework for dealing with 
dimensionality effects in ideal quantum gases, such as anomalous physics in d = [3] 
and a remarkable thermodynamic equivalence (discovered by Lee) between the Bose- 
Einstein and Fermi-Dirac gases in d = 2 |5,7J. There are suggestions that negative 
dimensionality, i.e. d < 0, might be of some theoretical interest [6j[8]. 

Motivated by the work of Lee [IHS] on the thermodynamics of quantum gases in 
terms of polylogarithms and, in particular, by his paper on polypseudologarithms, 
i.e. Li s (z) for s = 0, —1, —2, . . . , etc. [51 Sec. 3], we have begun a systematic study 
of Li s (z) of non-positive integral order s. Various known results regarding these 
functions, mainly obtained in different contexts and currently scattered throughout 
the literature, have been brought together here along with many new results and 
insights and they all have been proved in a simple and unified manner. However, 
it is neither intended nor attempted to provide a detailed historical account of this 
topic, nor to attribute the results to the original authors. Much attention has been 
given to considering all possible ways of computing Li s (z) for s = 0, — 1, —2, . . . , and, 
in particular, to the problem raised by Lee [6j Sec. 3] concerning the existence of a 
general closed-form expression for these functions. It has turned out that there are, 
in the mathematical literature, several such expressions involving either the Stirling 
numbers of the second kind or the Eulerian numbers. In Sec. 4 we shall deduce a 
new explicit closed-form expression involving the Carlitz-Scoville tangent numbers. 



The polylogarithm functions (polylogarithms, or polylogs for short) of order s 
and argument z, Li s (z), are defined by 



whenever this integral converges, i.e., 5?(s) > 0, ?R.(z) < 1, and elsewhere by analytic 
continuation. It is assumed that log(l/t) has its principal value and, evidently, there 
is branch cut from z = 1 to oo. In the important case where the parameter s is an 
integer, Li s (z) will be denoted by Li n (z) (or Li_ n (z) when the parameter is negative). 
Li n (2;), n S N := {1, 2, . . . , }, are the (classical) polylogarithms of order n (i.e. the 
nth order polylogarithms). The special case n = 1 is the ordinary logarithm (or 
monologarithm) Lii(z) = — log(l — z), while the cases n = 2,3,4, ... , are known, 
respectively, as dilogarithm, trilogarithm, quadrilogarithm, etc. 



2. POLYPSEUDOLOGS, THEIR PROPERTIES AND GENERATION 
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For more details and an extensive list of references in which polylogarithms appear 
in physical and mathematical problems we refer the reader to Maximon [9j. The 
nth order polylogarithms are thoroughly covered in Lewin's standard text [10], while 
many formulae involving Li s (z) can be found in Erdelyi et al. [114 pp. 30-31] and 
Prudnikov et al. p~2j pp. 762-763]. For some new results and applications see [i~3rfi~5] . 

Following Lee [6], the functions Li_ n (z), n G No, are here referred to as the 
polypseudologarithms (or polypseudologs) of order n and we first summarize their 
main properties and later thoroughly consider different ways of their computation. 
Polypseudologs exhibit several very important properties which are common to all 
polylogarithms. (i) The recurrence relation for poly logs can be rewritten as follows 

z^U^z) = U_ {n+1) (z) (neN ). (2.1) 
(ii) Polypseudologs satisfy the following particularly simple inversion relation 

Li_ n Q = (-l) n+1 Li_ n (z) (n € N). (2.2) 

(in) Duplication relation (or quadratic transformation) 

Li_ n (z) + Li_ n (-z) = 2 1+n Li_ n (z 2 ) (2.3) 

is satisfied, (iv) If \z\ < 1, polypseudologs have the series expansions given by 
X^feLi k n z n . In addition, hi_ n (z) have some additional, to these functions, specific 
properties, (v) Lio(z), Li_i(#), Li_2(^) 3 • • • , are rational functions, (vi) They have 
a pole of order n + 1 at z = 1. (vii) The values of polypseudologs at z = — 1 are 
related to the values of the Riemann zeta function at negative integers and are 
expressible in terms of the Bernoulli numbers B n 

Li_ n (-1) = (2 1+ ™ - 1) C(-n) = (1 - 2 1+ ") (n € N). (2.4) 

n + 1 

Hence, Li_ 2n (-1) = and Li!_ 2n (-1) = (1 - 2 2n )B 2n /(2n) since B 2n+ i = 0. (urn) 
Li_ n (z) are, for all orders, factorable by 2; and also by (z + 1) if n is an even number. 

Lee [HI Sec. Ill], by applying the recurrence relation to Lii(z) = — log(l — z) and 
repeating it over and again, obtained the Li_ n (z) to order n = 8: 

Li (z) = z/(l - z), 
Li_i(z) = z/(l - zf 
Li_ 2 (z) = z(l + z)/(l-zf, 
Li_ 3 (z) = z(l + 4z + z 2 )/(l-z) 4 , 
Li_ 4 (z) = 2(1 + z)(l + lOz + z 2 )/(l - z) 5 , 

Li_ 5 (z) = z(l + 26z + 66z 2 + 26z 3 + z 4 )/(l - zf, 

Li_ 6 (z) = z(l + z)(l + 56z + 246z 2 + 56z 3 + z 4 )/(l - z)\ 

Li_ 7 (z) = z(l + 120z + 1191z 2 + 2416z 3 + 1191z 4 + 120z 5 + z 6 )/(l - zf, 

Li_ 8 (z) = z{\ + z){l + 246z + 4047z 2 + 11572z 3 + 4047z 4 + 246z 5 + z 6 )/(l - zf. 

Clearly, in this way one can find the poly log of any desired lower order. In this 
section we consider a half-dozen additional ways for obtaining polypseudologs and 
majority of them are the closed-form formulae. 
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The second method to generate Li_ n (z) will next be described. Note that Li_i (z) = 
( z ~ttz) YlT=o zk i where we first differentiate and then multiply by z, and, if we apply 
the operator z4- n times, which we denote by the symbol fzj-) , to the summation 
formula Yl 



k=0 



[1 - z)- 1 , we have pH p. 364, Eq. (2) 



LL 



dz 



(n G N). 



(2.5) 



Moreover, it will be shown below that this relation leads to a closed- form formula 
involving the Stirling numbers of the second kind. 

An alternative way of generating hi_ n (z) for any n would be to make use of the 
generating function method, i.e. to generate {Li_ n (z)}^° =1 from a single function 
of two variables G(z, t) by repeated differentiation of that function. It is fortunate 
that there are several such functions (of which (2.6a) and (2.6b) could be found in 
the literature; see [2Ql p. 987] and [HJ P- 152, Eq. (19)]) 



°° t n 1 

V — Li_ n (z) = 

L — ' n\ 1 — ze 1 



n=l 

oo 



t" 



-nU 



n=0 

oo 



e L — z 



E^(-ir 1 LMn-i)W=log(e* 



n=2 

oo 



^ t n T . . . 11 + ze* 

> — U-n(z) = - 

n\ z 1 - ze 1 

n =1 

so that, for n G N, we have 



Li-nO) 



d^ 
dt n 



1 



z e' 



(-l) n z 



t=0 



d^ 



1 



171+1 

(-l) n T-ZT log(e* 



t=0 



2 dt n 



z / t=o 
1 + ze 1 



1 



z e 



t=o 



(2.6a) 
(2.6b) 
(2.6c) 
(2.6d) 



(2.7) 



Observe that it is in fact enough only to demonstrate that (2.6a) is the generating 
function for Li_ n (z); the function in (2.6b) follows from (2.6a) (as well as vice versa) 
by appealing to the inversion property in (2.2), while, the functions given by (2.6b) 
and (2.6c) could be transformed into each other by the recurrence relation (2.1). It 
is interesting that both functions 1/(1 — ze*) and (1/2)(1 + ze*)/(l — ze 1 ) have the 
same first derivative e*/(l — ze t ) 2 , thus we have that (2.6a) and (2.6d) are equivalent. 
Regarding (2.6a), let ze 1 < 1, where \z\ < 1, then, since Sfc°=o( 2:et ) fc = V(l — ze *); 
it is obvious that nth derivative of 1/(1 — ze 1 ) with respect to t evaluated at t = 
gives lA- n (z). 

Another way to compute Li_ n (z) stems from the fact that all the derivatives 
indicated in (2.7) are readily expressible analytically in a closed form. Indeed, if we 
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use the known result ( [2D p. 152] and [221 p. 32]) 



t=o v ' ' dt k 

k=l 



t=l 



where S(n, k) stands for the Stirling numbers of the second kind defined by means 
QUI P- 204, Eq. (lb)] 

k) = ± £(-D fc - (f)r = i £(-1)' (J) (* - *)", (2.9) 

' i=Q V 7 ' £=0 V 7 

(also see (2.12) below) and $(i) is any infinitely differentiable function, then, start- 
ing from the corresponding expressions in (2.7), we arrive at the following explicit 
formulae valid for n G N (for (2.10b), see [211 P- 152]) 

^W = tf^, (2-lOa) 



fe=l 



Li _ nW = £ HgW . ( , 10c) 

fc=0 v ; 

We remark that (2.10a) is obtained by (2.7a) as well as by (2.7d), and, furthermore, 
this formula also follows from (2.5) upon using [171 p. 364, Eq. (3)] 



m(z) = Y J Z k S{n,k)—±± (n € N), (2.11) 

/ / 1 



d 

k=l 



^f(z) being any infinitely differentiable function. It should be noted that it is easy 
to prove the formulae (2.8) and (2.11) by induction on n. For n = 1, they are 
obviously true and all that is needed for the induction step is the recurrence relation 
for S(n, k) pH p. 208, Eq. (3a)] 

S(n, k) = k S(n — l,k) + S(n — l,k — 1), 

S(n, 0) = 5(0, k) = and 5(0,0) = 1 (k, n G N). (2.12) 

Further, in view of the value of S(n,k) explicitly given by (2.9), it is straight- 
forward, starting from (2.10), to obtain another group of the closed- form formulae, 
this time involving double finite sums (for (2.13c), see [23 p. 988, Eq. (7)]) 



fc=l£=0 
n k 



J J (1-z) 

k=l£=0 \ / V I 

k=oe=o v y y ' 
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Next, it is not difficult to show that, in terms of the Eulerian polynomials and 
numbers, A n (z) and A(n, k), we have [Ml p. 245, Eq. (5n)] 

LM*) = (T^TT = (T^TT pMn,k)z k . (2.14) 

The Eulerian polynomials, A n (x), are usually defined by the generating function [161 
p. 244, Eq. (5i)] 

°° + n 1- x 

^^'r i-^i (2 - 15) 

n=0 

and they are polynomials in x of degree n, explicitly given by 

n 

A n (x) = J2 A (n,k)x k , (2.16) 

fc=0 

where the coefficients A(n, k) are the Eulerian numbers \16\ p. 242, Eq. (5f )] 

A(n,k) = Y J (-iy( n ^ 1 )(k-() n . (2.17) 

e=o \ * / 

Observe that the Eulerian polynomials and numbers are not to be confused with the 
Euler polynomials E n (x) and numbers E n . 

To prove (2.14) we may proceed as follows: upon replacing t by tj (1 — x) in (2.15) 
and dividing both sides of the obtained expression by (1 — x), the right-hand side of 
(2.15) then becomes 1/(1 — ace ) which is the same as the right-hand side of (2.6a), 
and hence (2.14). 

However, there is one more formula: by (2.14) and (2.17) we get ( [181 p. 325], |19] 
and [201 P- 987, Eq. (3)]) 

= tttVh t B-tf ( n r) ( * - i)n zK (2 - i8) 

V i fc=0£=0 V J 

3. New explicit closed-form formula for polypseudologs 

In this section we shall deduce the following new explicit closed-form formula for 
polypseudologs Li_ ra (z), n € N, 

(-i) L«/-u T(n , d + ^ LiliL r(n + !, fc) r i±£ ^ 

(3.1) 

where T(n, fe) stands for the tangent numbers (of order k) or the higher order tangent 
numbers defined by |23j p. 428] 

tan fe (t) = T{n, k) — (k € N). (3.2) 

n=fc 

Our derivation of (3.1) goes through two steps. The first step is to show that 
Li_ n (z) could be expressed in terms of the derivative polynomials for tangent, P n (x), 



Li-n(s) 



2 n+l 
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which are given by means of the exponential generating function 

P(Xj t ) := X + tanW , = Y P n (x) - t (3.3) 
v ' ; 1 -xtan(t) ^ K 1 n\ x ' 

v ' n=0 

or, equivalently, by the recurrence relation [24] 

P {x) = x, P n (x) = (l + x 2 )P' n _ 1 {x) (neN). (3.3*) 

The second step is to make use of a new formula for P n {x) 

n+l . 

P n («) =T(n,l) + ^-T(n + l,A;)^ (neN), (3.4) 
fe=i 

which shall be proved separately at the end of the section. 

We are now ready to proceed with the derivation of (3.1). First, consider the 
generating function for the derivative polynomials P n (x) in (3.3). Upon multiplying 
both sides by i/2, i := replacing x by — i(l + z)/{l — z), replacing t by it/2 

and noticing that tan(zt/2) = — i (1 — e*)/(l + e*), we obtain, after simple algebra, 
the following 

n =1 

Now, if we compare this with the generating function for Li_ n (z) given by (2.6d), 
then, clearly, we have 

/ 7 \ n+l / 1 i _i_ z \ 
Li_ n (z) = (-j P n (-.-L±) ( nG N). (3.6) 

Secondly, after combining (3.6) and (3.4), we get a formula similar to the required 
formula (3.1) and what remains is to evaluate the powers i n+1 and % n ~ k+1 which 
respectively appear as multiplication factors before T(n, 1) and T(n + 1, k). At this 
point it is helpful to examine the tangent numbers in more detail (see Table 1 and 
Eq. (3.2)) and notice that T(n,k) ^ only when 1 < k < n and either both n 
and k are even or both n and k are odd. In other words, T(2m, 2r + 1) = and 
T(2m + 1, 2r) = (m, r £ N ). 



Table 1. Tangent numbers T(n,k) 



n\k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


1 


1 


















2 





2 
















3 


2 





6 














4 





16 





24 












5 


16 





120 





120 










6 





272 





960 





720 








7 


272 





3 696 





8 400 





5 040 






8 





7 936 





48384 





80640 





40 320 




9 


7396 





168 960 





645120 





846 720 





362 880 
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This means that the terms involving T(n, 1) when n is odd and T(n + 1, k) when 
n — k is odd are the only surviving terms in the obtained formula. Hence, we only 
need to evaluate the powers z n+1 when n is odd and % n ~ k+1 when n — k is odd. This 
is how we have found the negative powers, so the proof of (3.1) is now complete. 

In order to prove formula (3.4) note that the generating function of P n {x) in (3.3) 
can be rewritten as 

oo oo 

P(x, t) = [x + tan(t)] i x tan(t)] fc = x + (l + x 2 ) Y x k ~ l tan fc (i) 

fc=0 k=l 

which, by making use of the definition of T(n, k) in (3.2) and the following elementary 
double series identities [25l p. 57, Eq. (2)] 

oo n oo oo oo oo 

Y Y a ( k ' n ) = Y Y a ( /c > n + k "> = YY °( fc > n )' 

n=lfe=l n=lfc=l n=lk=n 

becomes 

t n 
nl 

k=l n=k 
oo / n \ 

=x+ j2(i+x 2 ) i Yn^^x^ 1 ) -. (3.7) 

Tt — 1 \ — 1 / 

On the other hand, by the definition of P n (x) in (3.3) in conjunction with the 
recurrence relation for P n {x) (3.3*), we have 



P(x, t) = x + (1 + x 2 ) Y Y T ^ n ' fc ) 



P(x,t)=P (x) + Y Pn(x)-=X+Y{ l + X 2 )P'n-l{x)-. (3.8) 

z — ' n! z — ' n! 

n=l n=l 

and thus comparing (3.7) with (3.8) clearly yields 

n n+1 ^ 

P^W = J^r^fcjx*- 1 , thus P n (x) = P n (0) + YvT(n + l,k)x k . (3.9) 



fc=l k=l 

Finally, the desired result (3.4) follows from (3.9) since, by the definitions (3.2) 
and (3.3), we have P n (0) = T(n, 1). 

4. Concluding remarks 

It is interesting to note that several special functions, such as the Riemann zeta 
function C( s )> the- Hurwitz zeta function C( s : a ) an d the classical gamma function 
r(z), when the argument of a function takes on non-positive integer values show 
an abrupt changes in properties which cannot be easily explained. For instance, at 
z = 0, —1, —2, . . . , the function T(z) becomes singular, while £(s) takes on rational 
values. On the other hand, £(s, a), for s = 0, —1, —2, . . . , reduces to a polynomial 
in a which is expressible in terms of the Bernoulli polynomial B n (x) as follows 

a-n,a) = -^±^ (n G No). 
n+1 

Similarly, the polypseudologs considered in this paper, i.e. the polylogs Li s (z) for 
s = 0, — 1, —2, . . . , are, in fact, rational functions in z. 



Polypseudologarithms revisited 



9 



Lee [1] already made use of the nil-log Liq(z) = z/(l+z) in his study of the statisti- 
cal mechanics of ideal particles in null dimension and discovered several anomalous 
physical phenomena, among them probably the most remarkable one is that the 
bosons are not confinable. It remains to be seen whether it would be possible, as 
it is expected, to employ the polypseudologs Li_i(z), Li_2(z), Li^(z), . . . , in some 
kind of extension of Lee's formalism to the case of negative physical dimensions. It 
is hoped that this paper will contribute to a better understanding of the analytical 
properties of these functions and stimulate further work in the field. 
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